Comment on "Time-reversal symmetry breaking superconductivity" 
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It is pointed out that erroneous Bardeen-Cooper-Schrieffer model equations have been used by 
Haranath Ghosh in his recent treatment of time-reversal symmetry-breaking superconductivity. 
Consequently, his numerical results are misleading, and his conclusions are not to the point. 
PACS number(s): 74.20.Fg, 74.62.-c, 74.25.Bt 
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Some recent studies provide increasing evidence that 
the pairing symmetry of some of the cuprates at low 
temperatures allows an order parameter in a mixed sym- 
metry state. At higher temperatures, below the critical 
temperature T c , the symmetry of the order parameter is 
of the d x 2_ y 2 type. At a lower temperature there could 
be an admixture of a minor component, such as d xy , on 
the predominant d x i_ y i symmetry. This general time- 
reversal symmetry breaking order parameter has the form 
d x 2_ y 2 + exp(i9)dxy, where 9 is the mixing angle. 

Recently, Ghosh pj| presented a theoretical study of su- 
perconductivity for this mixed-symmetry case based on 
the Bardeen-Cooper-Schrieffer (BCS) equation. In this 
comment we point out that the coupled equations used 
by him for the two components of the order parameter are 
erroneous and present a rederivation of the appropriate 
equations. This comment also applies to latter investiga- 
tions where Ghosh further used the erroneous equations 
in subsequent studies of (1) superconductors of mixed or- 
der parameter symmetry in a Zeeman magnetic field 
and (2) pairing symmetry and long-range pair potential 
in a weak-coupling theory of superconductivity |3|]. 

We use the two-dimensional tight binding model as in 
Ref. jlj] . The effective interaction after including the two 
appropriate basis functions is taken as 



"Vl?7lk?7lq - V2?72k?72q, 



(0.1) 



such as, d x 2_ y 2, s (r/ = 1), s x 2 +y 2 (r/q = cosq x + cosq y ), 
and s xy (jjq = cos q x cos q y ), and the present discussion 
equally applies to mixtures involving such orthogonal 
states. On the continuum the s-wave angular function 
Ci(^) — 1j an d the d-wave circular harmonics ^(0) = 
cos(2</>) and ^((j)) = sin(20) satisfy the trivial orthogo- 
nality relation 



2tt 



(0.3) 



One passes from the lattice to continuum description by 
replacing (a) the sum over q by an integral over </> and 
(b) the functions rjq by the circular harmonics £(4>). 

Although Ghosh jlj] considered the BCS model at a 
a finite temperature, we consider its zero temperature 
version, which is enough for our purpose: 



(0.4) 



with Eq = [(e c 



M) 2 



jAqj 2 ] 1 / 2 , where e q is the single- 



particle energy and fi is the chemical potential. The order 
parameter has the following general anisotropic form: 



Al77l q + CA 2 772q, 



(0.5) 



cosq y corresponds to d x 2_ y 2 sym- 

smq y corresponds to d xy symme- numbers. ' If we substitute Eqs. (|0 



where r/iq = cos q x - 
metry, 77 2q = sin^ 

try. The orthogonal functions ?7iq and ?72q are associated 
with a one-dimensional irreducible representation of the 
point group of square lattice C^ v Q and can be consid- 
ered appropriate generalizations of the circular harmon- 
ics cos(2</>) and sin(20) incorporating the proper lattice 
symmetry. The orthogonality condition of these func- 
tions is 



where C = exp(i0) — (a + ib) is a complex number of unit 
modulus | C | 2 = 1 and a = cosO and b = smO are real 



and ( |0.5| ) into the 

BCS equation ( p.4[ ), for orthogonal functions r/iq and ?72q, 
one can separate the resultant equation into the following 
components 



VlqV2q = 0. 



(0.2) 



q 

(a + ib)A 2 = V 2 J2 



»7lq[AlJ7lq 



i6)A 2 r?2q] 



2Eq 

??2q[Al?7lq + (fl + ib)A 2 1]2q} 

2E a 



(0.6) 



(0.7) 



This orthogonality relation is readily verified as under the 
transformation q x — » — q x (or q y — > — q v ) rj2q changes sign 
and ?7iq remains unchanged. Although, the proof could 
be slightly different, similar orthogonality relation exists 
between the basis functions of states on square lattice, 



Equations ( |0.6| ) and (|0.7| ) have solution for real Ai and 
A2, when the complex parameter C is either purely real 
or purely ima gina ry. 

Equations (CK6) and (CV7) can be substantially simpli- 



fied for a purely imaginary C, e.g., for C — i or a = 



1 



and 6=1 



7r/2). In this case for real compon ent s 



Ai and A2, the real and imaginary parts of Eqs. 
and ( |0.7| ) become, respectively, 



d__j) 



A 2 = V 2 



q 



?7 2 2 q A 2 
2£„ 



Here we have used the identity 



E 



?7lq?72q 

2E n 



= 0, 



(0.8) 



(0.9) 



(0.10) 



Jq = [(e, - ^) 2 + A 2 r; 2 q 



which holds in this case as i? g 

A^lq] 1 / 2 is invariant under transformation q x — » — (fo 
or under — > — (fo, whereas under either of these trans- 
formations 772 q changes sign and ?7iq remains unchanged. 
Henc e usin g the integrat ion of q, o ne c an establish iden- 
tity (3.10). Equations (0.8) and (3.9) have been used 
in the study of the mixed-symmetry states of types 
d x 2_ v 2 + id xy and d x 2_,p + is H. 



Equations (3.6) and (3/7) also lead to a simple form 
for a purely real C, e.g., for C = 1 or for a = 1 and 6 = 
(8 = 0). However, in this case E q = [(e q — fi) 2 + (Airyiq + 
A 2 ?72q) 2 ] 1 ^ 2 contains cross terms of the type ryi q ?72q, and 
is not invariant under transformati on q x — > — q x or under 
9y — * — #y Consequently, Eq. ( p. 10 ) is not satisfied 
and coupled angular terms will be present in the BCS 
equa tion. In this case for real Ai and A 2 , Eqs. ( |3~6| ) and 
(3.7) become the following set of coupled equations 



Tyi q [Ai?yiq + A 2 r? 2q ] 
2E n 



(0.11) 



A2 = ^w^m (012) 



In t he ca se of a general mixture, e.g., a ^ and 6 ^ 0, 
Eq. (3.10) is not valid. How ever, f or r eal Ai and A 2 , 
one can break up Eqs. (D.6) and ( |0.7j ) into their real 
and imaginary parts, e.g., into four cou pled equations 

As Eq. 



for two unknowns Ai and A 2 . As Eq. ( |0.10 ) does not 
hold in this general case, the four coupled equations are 
consistent only if A 2 = 0, or Ai = 0, which corresponds 
to no coupling between the two components. Hence the 
permissible values for the mixing angle 9 are 0, 7r/2, 7t, 
and 3n/2. 

In his study Ghosh [Q implicitly assumed Eq. ( O.lCj ) 
to be valid in all the cases discussed above including (i) 
the case of a general C with a ^ and 6 ^ 0, and (ii) 
the case with a = 1 and 6 = 0. Consequently, he arrived 
at the wrong equations ( |0.8D and ( |0.9| ) for a general C, 
which he used in his numerical treatment, specifically for 
mixing angles 9 = 0, and 7r/4 (see Eq. (9) of Ref. (p, 



Eq. (6) of Ref. §, and Eq. (5) of Ref. §). For 9 = 
he used inappropriate equations and for 9 = it/ 4 there 
should not be any mixing. 

It is in tere sting to recall that using the orthogonality 
relation (D.3) on the continuum Musaelian et al M de- 
rived the BCS equations for the mixed-symmetry states 
s + d and s + id. In agreement with the present comment 
and in contradiction with the investigation by Ghosh 
|l|-^| they (a) confirmed the existence of mixed-symmetry 
states for the mixing angles 9 = and ir/2 only and (b) 
reported the BCS equation for the s + d state, which is 



structurally quite similar to Eqs. (0.11) and (0.12) above. 
As the study of Musaelian et al || referred to the contin- 
uum, in that work the discrete sum over q was replaced 
by the integral over <f) and the functions 77 replaced by 
the circular harmonics £. 

Rec ently, we used the correct equations (|0. 11 ) and 
( |0. 12 ) for a description of the d x 2_ y 2 + d xy symmetry 
case [0, which corresponds to 9 = above. The quali- 
tative feature of the temperature dependence of the A's 
in that study is quite distinct from the erroneous results 
obtained by Ghosh |l],§ by using the inappropriate equa- 
tions. 
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